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An intermediate mass black hole (IMBH) may have a dark matter (DM) mini-halo around it 
. and develop a spiky structure within less than a parsec from the IMBH. When a stellar mass 

(N . object is captured by the mini-halo, it eventually infalls into such an IMBH due to gravitational 

I wave back reaction which in turn could be observed directly by future space-borne gravitational 

. wave experiments such as eLISA/NGO. In this paper, we show that the GW detectability strongly 

depends on the radial profile of the DM distribution. So if the GW is detected, the power index, that 
' is, the DM density distribution would be determined very accurately. The DM density distribution 

^S) ■ obtained would make it clear how the IMBH has evolved from a seed BH and whether the IMBH 

has experienced major mergers in the past. Unlike the gamma ray observations of DM annihilation, 
GW is just sensitive to the radial profile of the DM distribution and even to non-interacting DM. 
Hence the effect we demonstrate here can be used as a new and powerful probe into DM properties. 



INTRODUCTION 



A large number of astrophysical and cosmological observations provide convincing evidence for the existence of dark 
matter (DM). The origin and nature of DM remain largely unknown, and are among the most challenging problems 
J — ■ in current cosmology and most likely in particle physics. 

^\ [ Recently, the distribution of DM around a black hole (BH) has been under discussion in the context of indirect 
• searches for DM annihilation signals with gamma ray observations. Gondolo and Silk [l[ first suggested that the 
adiabatic growth of a BH creates a high density DM region, called the "spike" , which enhances the DM annihilation 
■ rate. Subsequent work showed that the existence of a DM spike around a supermassive black hole (SMBH) turns 
[ out to be unlikely when one considers the effects of major merger events of the hostgalaxies [ll, off-center formation 
of the seed BH and scattering of dark matter particles by surrounding stars On the other hand, a DM 

"mini-spike" around an intermediate mass black hole (IMBH), with a mass range between 10^ to lO^M©, may survive 
if the IMBH never experienced any major mergers 0,0]; as is expected to be the case for the many IMBH that have 
failed to merge into a supermassive BH. 
d • The existence of such a spike structure is strongly dependent on the details of BH formation and the history of major 
mergers, which are far from clear. In this Letter, we propose that future gravitational wave (GW) experiments can 
be used to probe the DM distribution around BHs. The existence of the dense DM region changes the gravitational 
potential and affects the orbit of an object around the BH. We consider GWs from the coalescence event of a compact 
binary consisting of a small mass object and an IMBH, and evaluate the modification of the GW signal by the existence 
of a DM spike associated with the IMBH. Such an event may be observed by future space-based interferometers such 
as eLISA/NGO [l] and DECIGO 0. We further discuss whether the eLISA/NGO experiment is sensitive to the 
modification of the signal by the DM spike. 

Note that, while gamma ray observations can find the signal of DM annihilation if DM is a weakly interacting 
massive particle, the observation of GWs is just sensitive to the gravitational potential of the DM halo and applicable 
for any type of DM. Therefore, future GW experiments offer a unique opportunity for testing the existence of the DM 
spike around BHs. 

Let us describe the radial profile of the DM spike by a single power-law p cx r~" assuming a spherically symmetric 
distribution of DM. The adiabatic growth of the BH produces a dense spike in the inner region of the mini-halo 
within a radius of Vgp ^ 0.2r/i, where is the radius of gravitational infiuence of the BH defined by Af(< Vh) = 
4tt J^'' plry^dr = 2Mbh, with A/bh being the BH mass [J]. The final density profile of the spike depends on the 
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power-low index a-mi of the inner region of the initial mini-halo as, a = (9 — 2aini)/ (4 — aini) [H, [13 • If '^^ assume the 
Navarro, Frenk and White (NFW) profile [13] for the initial condition (ajn; = 1), we get a = 7/3. Very steep slope is 
generically predicted as we find 2.25 < a < 2.5 for < aini < 2. 

In summary, in this paper, we assume the DM distribution of a mini-spike is described by 

p{x) = psp {-^"^ {rmin <r < r^p), (1) 

where psp is the normalization of the DM density. For an IMBH with the mass of Mbh = 10'^ M© and the total mass 
of the DM mini-halo of Mhaio = lO^M©, we get psp = 379 Mq/pc^ and Vgp = 0.33pc. Beyond the spike radius Vgp, 
the DM distribution obeys the NFW profile with the concentration parameter c = 6.6 estimated based on the fitting 
formula given, e.g., in The minimum radial distance is taken to be rmin = nsco where nsco is the Innermost 
Stable Circular Orbit (ISCO) given by rmin = nsco = GGMbh/c^ ■ Throughout the paper, we restrict a to below 3, 
since Vgp is nearly the same as rmin and the mini-spike virtually vanishes when a > 3. 



FORMULATION 
GWs from binary inspiral 

Let us consider gravitational waves from a binary system consisting of an IMBH with a mass of Mbh ~ IO'^Mq and 
a compact object with a mass of IMq. For simplicity, we make the following idealization. First, we treat the star 
as a test particle and we call it a "particle" in the following. Second, we assume that the DM density is unperturbed 
even when the star orbits in the DM halo. Gravitational heating of the DM mini-spike due to the particle may be 
noticeable about within the Hill sphere of the particle because of the gravity of the central IMBH. In the case of our 
IMq-10'^Mq binary, the Hill radius is 10% of the orbital radius and we ignore possible heating effects in the first 
order approximation. Then, the equation of motion for the particle is written as 

d^r _ GM.a F P 
where I is the angular momentum of the particle per its mass, and M^s and F is 

Meff ={ 3-a , F={ 3~a (3) 

Q yT <^ '"min j 

In the first term of the right hand side of Eq. the DM mini- halo modifies the effective mass of the central IMBH. 
The second term contains information of the DM halo radial distribution. The third term represents a centrifugal 
force. If we assume that the second term is much smaller than the first term. 
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we can treat the term which involves information on the DM halo as a perturbation, and expand equations in powers 
of e, which is a dimensionless parameter depending the power index a. 

When the particle stably orbits around the IMBH at a constant radius R, the left hand side of the equation of 
motion vanishes. In this case, the GW waveforms are given by 



1 iGuuj^R'^ 1 + cos'^ L . 
h+ = ^ ^ cos 2c.,i 4 



Aix = 7 cos i sm (2a;st) (5) 



to the lowest order approximation where l is the inclination which is the angle between the normal to the orbit and 
the line-of-sight, and 2los is the GW frequency. 
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Waveforms including GW back-reaction 

Next, we include the effect of tlie GW back-reaction within the hnearized theory of Einstein's general relativity. The 
orbital radius and frequency is no longer constant, because GW radiation energy i?GW is taken from the rotational 
energy ii^orbit of the particle. The relation between the orbital radius R and the time t is given by (e.g., [ist). 



orbit 



dt 



dEQ\fi/ 
di~' 



(6) 



Using this relation, we can compute the orbital frequency ujg and i? as a function of time. To include the GW back- 
reaction in the GW waveforms, we replace the constant parameters and R in Eqs. ^ and ([5]) by time-dependent 
functions uJs{t) and R{t). 

The stationary phase approximation then enables us to obtain the GW waveforms in Fourier space expanded in e, 
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where tc is the value of retarded time at coalescence, $o is the value of the phase at coalescence, Mc = /i^/'^M^g^ is 
the chirp mass, and ^' is the phase of the GW waveform. These expansions are valid for the frequency / for which 
higher order terms are negligible. 

In Eq. ([9]), the phase of the GW is modified by the presence of the DM, which is expressed in powers of e. Since 
GW interferometers are very sensitive to the phase of the signal, this phase difference is crucial for distinguishing the 
existence of the DM spike. In Fig. [l] we plot the phase difference caused by the DM halo, taking into account 
terms up to second order in e. We see A^ increases for low frequencies and for large a. This can be explained by 
the fact that the orbit of the object is affected only by the DM mass inside the orbital radius. More phase difference 
is produced when the inner mass is large. As shown in Fig. [2l the enclosed DM mass increases as the radius or a 
increases. Since a low frequency of the GW corresponds to a large orbital radius, a large phase difference is produced 
at low frequencies. Large values of a means a steep density distribution and leads to a large inner mass, which also 
results in a large phase difference. 



OBSERVATION OF GWS 



Matched filtering 



Let us discuss if this effect is testable by future GW experiments. The search for GW signals is performed by 
matched filtering analysis, in which one correlates detector output with theoretical template. The signal-to-noise 
ratio is defined by 



^1 =4^ 

N 



df 



hif)h*Af) 
S{f) 



df 



/ini 



ht if) 



(11) 



4 





radius r (pc) 



FIG. 2: Mass of DM halo within orbital radius r. Steeper density distributions contain more DM mass within the radius r. 
The solid line is a = 1.5, the dashed line is a = 2.0 and the dot-dashed line is a = 7/3. 



where h{f) is the GW signal coming to the detector, ht{f) is the template, S{f) is the spectral density of the detector 
noise, and fi^i is the frequency of the inspiral GW when the observation started. In the following example, we assume 
the eLISA experiment, whose noise spectrum is given in Ref. @. 

In Eq. the numerator is a noise- weighted correlation between the template and the true signal, and the 

denominator is the renormalization factor. When the template matches the true waveform, S/N is maximized. Thus, 
S/N is an indicator to tell us whether the waveform of the template is present in the detector or not. 



Detectability of the effect of a DM halo around an IMBH 



Let us consider an observation of GWs from the IMq particle inspiralling into the 10"^ Mq IMBH, which would be 
detectable by the eLISA experiment. We assume this binary is surrounded by a DM mini-spike whose distribution 
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FIG. 3: P/Q against power index a. Three different curves show P/Q for three different values of initial frequency /ini, namely 
differential observation time. The solid line is for /ini = 0.1 Hz, the dashed line is for /ini = 0.01 Hz and the dot-dashed line is 
for /ini = 0.001 Hz. 



is given by Eq. ([T]). In this setup, a frequency integration from /ini = O.OlHz corresponds to a 4.5 year observation 
until the coalescence (corresponding roughly to the expected eLISA frequency band and observation time) . Note that 
when /ini = O.OlHz, the particle is about 10^^ pc away from the IMBH and well within the mini-spike. 

In Fig. [3j we show how much the S/N is degraded when one applies a template predicted without considering 
the DM effect on the signal with the DM effect. The vertical axis represents a degradation rate P/Q, where P is 
S/N calculated assuming a template of a waveform without the DM effect (e — 0) and Q is S/N calculated with a 
template including the DM effect up to the second order. If the effect of the DM is small, there is little difference 
between the two templates and P/Q becomes 1. Conversely, if DM potential induces significant phase difference, the 
value of P decreases, since the template and the signal have less correlation. 

As discussed in the previous section, the phase difference becomes significant for large a, and, from Fig. [3l wc find 
P/Q largely deviates from 1 for a > 2. This result indicates that GW observation can distinguish whether a DM 
halo of a > 2 exists around the IMBH. In order to extract inspiral signals under the effect of a DM halo, we must 
prepare templates including the DM effect. In Fig. [3l we also plot the cases for different initial integration frequencies, 
which corresponds to different observation time. Since the phase difference becomes larger at low frequencies, P/Q is 
suppressed for smaller value of a when one observes a longer time period. The peak seen at a = 2.5 originates from 
the zero crossing of the first term of Eq. ([T0| . 



CONCLUSION 



We have demonstrated a method to probe the DM distribution around an IMBH by using GW direct detection 
experiments. Considering a GW inspiral signal from inspiral of a compact object falling in an IMBH, we have 
computed how the GW waveform is modified by the gravitational potential of the DM halo. Thanks to the fact that 
a GW interferometer measures the phase of the signal with very good accuracy, we found that a GW experiment such 
as eLISA/NGO is sensitive to the phase shift caused by the DM potential. Accordingly, GW observation may provide 
a new observational probe for the DM distribution near BHs and would be helpful for testing the existence of a dense 
DM minihalo, the so-called DM spike. This may even offer hints to the formation history of BHs, since formation of 
DM spikes strongly depends on how BHs evolved. 

In future work, we plan to extend the investigation for different values of mass distribution parameters, such as the 
power index a and the mass of the compact objects and the halo. Wc will also estimate to what degree future GW 
experiments can determine the mass distribution by computing expected errors on the parameters. 

Finally we remark on the dynamical friction effect. In this paper, we have ignored dynamical friction. Even if GW 
back reaction is ignored, the dynamical friction alone makes the smaller compact object infall into the IMBH. Using 
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the Chandrasekhar formula (e.g., [14|), we found however that this effect produces negUgible phase difference when 
computing the SN degradation rate. 
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